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1. Introduction 

Dynamics near the critical point can be classified depending on gross variables. 
For example, if the critical dynamics is described by a conserved order parameter, 
the dynamical universality class belongs to model B in the classification scheme of 
Hohenberg and HalperinpQ. If it is described by a one-component conserved order 
parameter and a three-component conserved quantity, the class belongs to model H. To 
see the dynamical universality class, it is important to know the gross variables that 
contribute to the critical dynamics. 

Recently, Son and Stephanov discussed the coupled dynamics near the QCD critical 
point phenomenologically|2]. They concluded that if we ignore the mode coupling to the 
energy density and the momentum density, the dynamical universality class of the QCD 
phase transition belongs to model B, and if we take them into account, it is changed 
to model H. However, there are possibilities that have not been discussed by them. 
1) They assumed that the time evolution of the gross variables (the chiral condensate 
and the net baryon number density) is proportional to thermodynamic forces defined 
by the derivatives of a Ginzburg-Landau free energy in terms of the gross variables. 
This approach will be adequate for deriving equations for non-conserved quantities, for 
example, the chiral order parameter. Then, we obtain the time-dependent Ginzburg- 
Landau (TDGL) equation. However, for conserved quantities like the net baryon number 
density, we should adopt the Cahn-Hilliard (CH) equation instead of the TDGL equation 
|3j- 2) To describe the critical dynamics of the QCD phase transition correctly, we should 
take into account the coupling of the chiral condensate and other gross variables, i.e., the 
net baryon number density, the energy density and so on. Then, it is known that there 
exist two channels that give rise to the mode coupling. One is the irreversible coupling 
where the mode coupling is induced through irreversible currents, and the other is the 
reversible coupling induced through reversible currents [TJ El IH El E|- The currents 
defined by thermodynamic forces are irreversible currents. Thus, Son and Stephanov 
discussed only the former possibility. However, it is known that the important mode 
couplings are usually caused by not the irreversible currents but the reversible currents, 
for instance, the liquid-gas phase transition, the superfluid transition of liquid helium 
and so on. 

In this paper, we have two purposes. First, we derive a phenomenological equation 
for the net baryon number density by using the CH equation. Then, the derive equation 
is different from that of [2,. Second, we show how the reversible current can affect the 
relaxation behavior of the chiral condensate. 

This paper is organized as follows. In Section 2, we derive the phenomenological 
equations of the chiral condensate and the net baryon number density. In Section 3, 
we discuss the irreversible coupling of the chiral condensate and the baryon number 
density. We derive the reversible coupling in the projection operator method in Section 

4. Then, we found that the chiral condensate is coupled to the chiral current. In Section 

5, we show that the reversible coupling can change the critical dynamics. Summary and 
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conclusions are given in Section 6. 



2. Phenomenological equations based on the TDGL and CH equations 

We adopt the same Ginzburg-Landau free energy / as used by Son and Stephanov0, 

F[a,n] = J d 3 x[^{d i0 -) 2 + bd^n + |(^n) 2 + V{a,n)\, (1) 

where a and n denote the fluctuations of the chiral condensate and the net baryon 
number density that is, a = qq — (qq) eq an d n = qrfq — {q , J q)eq, respectively. The 
potential term V(o~, n) is given by 

A C 

V(o~, n) = —a 2 + Ban + —n 2 + higher order terms. (2) 

The parameters A, B, and C satisfy the relation AC = B 2 at the critical point. 

The thermodynamic forces are defined in the following way. Usually, for this 
purpose, we should derive the Gibbs-Duhem relation including the chiral condensate 
and calculate the entropy production. On the other hand, the thermodynamic forces 
are assumed to be defined by the derivatives of the free energy in terms of a and n 
in [2] . In this paper, we assume that the phenomenological equations of the chiral 
condensate and the net baryon number density are given by the TDGL equation and 
the CH equation, respectively. From this assumption, we define the thermodynamic 
forces. 

First of all, we ignore the mixing between the chiral condensate and the net 
baryon number density. The chiral condensate is a non-conserved quantity. Then, the 
phenomenological equation of the chiral condensate is given by the the TDGL equation 
with noise, 

5F 

<i(x) = - J CT (x) + £ CT (x) = -Iss^ + &(x), (3) 

da 

where 7 SS denotes the Onsager coefficient associated with a(x). The last term £o-(x) 
represents the noise term of the chiral condensate. We do not give the concrete 
expression of the noise term because the following discussions are not affected by the 
noise. On the other hand, the net baryon number density is a conserved quantity. Then, 
the equation of the net baryon number is given by the CH equation with noise, 

8F 

n(x) = - VJ n (x) + f n (x) = -7n„V 2 — + £n(x), (4) 

on 

where 7„„ denotes the Onsager coefficient associated with n(x) and £ n (x) represents the 
noise term. 

Then, the currents of the chiral condensate and the net baryon number density are 
defined by 

5F 

J CT (x) = IssJ^, (5) 

5F 

J n (x) = 7rmV— , (6) 
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respectively. These currents are irreversible because the time reversal symmetry is 
broken. In nonequilibrium thermodynamics, irreversible currents are proportional to 
thermodynamic forces. Thus, the two thermodynamic forces are introduced by 
SF 

X a = °— = (A- aV> + (B- bV 2 )n, (7) 
SF 

X n = V— = V{(B-bV 2 )o- + (C-cV 2 )n}. (8) 
on 

It should be noted that the thermodynamic force X n defined here is vector, although 
the scalar thermodynamic force is introduced in j2]. 

Next, we will discuss the mixing of the chiral condensate and the net baryon 
number density. The general irreversible current is given by the linear combination 
of the thermodynamic forces. By using the two thermodynamic forces defined above, 
thus, the irreversible currents including the mixing of the chiral condensate and the net 
baryon number density are given by 

J CT (x) = 7 SS X CT + 7 sn X n , (9) 
J„(x) = 7„„X n + 7 ns X CT . (10) 

The irreversible current of the chiral condensate Q is scalar and hence the Onsager 
coefficient 7 sra must be vector. Similarly, the irreversible current is vector and hence 
the Onsager coefficient 7„ s must be vector. Then, by using the Onsager's reciprocal 
theorem, we assume 

Isn = Ins = aV '■ (11) 
Then, the off-diagonal Onsager coefficients have the q-dependences as is discussed in 

m 

However, in this paper, we adopt, what we call, the Curie principle (symmetry 
principle), where, for instance, it is forbidden that vector currents are induced by scalar 
thermodynamic forces [3 IHl E]- For instance, the thermal conduction current is not 
induced by chemical reactions because of the Curie principle. In our case, the irreversible 
current of the net baryon number density is not induced by X a . Accordingly, because 
of the Onsager's reciprocal theorem, the off-diagonal Onsager coefficients ^ ns and 7 sn 
vanish. Finally, the phenomenological equations are given by 

a(x) = - lss [{A - aV>(x) + (B - 6V 2 )n(x)] + &(x), (12) 
ra(x) = - 7n„V 2 [(£? - 6V 2 )a(x) + (C - cV 2 )n(x)] + £ n (x). (13) 

Because of the definition of the Ginzburg-Landau free energy, the equations still contain 
the cross terms that give rise to the mode coupling between the chiral condensate and 
the net baryon number density. 

It should be noted that if we ignore the coupling between a and n, it is easy to solve 
the linear equation and we can see that the relaxation time of the chiral condensate at 
the vanishing momentum is given by (7 ss yl) _1 . This result is consistent with the result 
of the van Hove theory [TU]. If the van Hove theory is exact, the dynamical critical 
exponent is completely decided by the temperature dependence of A, and hence takes 
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the same value as the static critical exponent. However, the experimental data shows 
the deviation from the static critical exponent. The deviation comes from the couplings 
between the order parameter and other gross variables, as we will see below. 



3. Irreversible coupling in the QCD phase transition 

In this section, we discuss the mode coupling induced by the irreversible current. 

After Fourier transformation, the Langevin equations in the limit of q — > are 
given by 

<Kq) = (-T«^ ~ ls S aq 2 )o-(q) + (-J SS B - -f ss bq 2 )n(q) + £ CT (q), (14) 
n(q) = j nn Bq 2 o-(q) + 'y nn Cq 2 n(q) + f„(q). (15) 

Here, we ignored the terms proportional to q 4 . 

Following the discussion of Son and Stephanov, we calculate the eigenmodes of the 
coupled equation. The eigenmodes are calculated by solving the equation 

-ls S A - j ss aq 2 + iu -J SS B - -f ss bq 2 

InnBq 2 InnCq 2 + iLO 



0. 



(16) 



Up to second order in q, we obtain the following two eigenmodes, 



. lnn (AC-B 2 ) 2 

Wi= -% q, (17) 



i '( r< \ — i -lnn{.AB — C j o 

^2= - ilssA + H7„„C - 7 ss a)g +i q. [II 



The corresponding eigenstates are 

\-A + 0(q 2 )J \ 0(q 2 

respectively. 

Afterward, Son and Stephanov developed the argument that the chiral condensate 
and the net baryon number density are not independent variables near the critical point 
and their linear combination gives one gross variable as is the case with the liquid-gas 
phase transition. 

However, if we take into account the possibility of the critical slowing down, the 
conclusion can be changed. From the results of the eigenmodes and eigenstates, one can 
easily write down the solution of Eqs. (fT4*|) and (fl3j) at low q, 

< cr(q,t) > = e-^ 2 V(q), (20) 

« n (q, t ) » = e-^< (n(q) - |a(q)) + |e-**V(q), (21) 

where <C S> denotes the expectation value for the noise with a suitable Gaussian weight 
function. It should be noted that there is a constraint for the initial state: n(0) = 0, 
because n(q) is the fluctuations. To represent dissipation, ^ SS A and 7 nn (y4_B — C 2 )/A 
must be positive. Furthermore, 7 SS A becomes smaller as the temperature approaches 
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the critical temperature because of the critical slowing down [TBI ITo] IT7| . Thus, it is not 
obvious whether we can conclude iu>2 is much larger than iui near the critical point. 
Then, the chiral condensate and the net baryon number density do not have a clear 
separation of time scales and cannot be reduced into one gross variable. 

4. Reversible couplings in the QCD phase transition 

In the previous section, we derive the irreversible currents of the chiral condensate and 
the net baryon number density in the framework of nonequilibrium thermodynamics. 
However, in general, the time-evolutions of hydrodynamic variables are dominated by 
not only irreversible currents but also reversible currents. The former yields dissipation 
and entropy is increased, and the latter yields macroscopic motion as in some cases of 
convection and entropy is not changed. This structure is general and does not depend 
on whether the variable is conserved or not. Thus, there can be other mode couplings 
caused by the reversible currents that have been ignored so far. As a matter of fact, 
it is known that the reversible coupling plays a dominant role in the critical dynamics: 
the most important mode couplings are ordinarily caused by not irreversible currents 
but reversible currents [IJ El 13 El El • In this section, we derive the reversible currents 
that can affect the critical dynamics. 

First, we consider the net baryon number density. The equation of continuity of 
the net baryon number density is 

ra(x) + VJ re (x) = 0, (22) 

where the net baryon number current is given by 

J re = _ ($ 7 V) e? . (23) 

Thus, the net baryon number density can be coupled to the net baryon number current. 
In short, the Langevin equation for the net baryon number density including the 
reversible current is given by 

n(x) = - VJ re (x) - 7nn VX n (x) + f n (x). (24) 

On the other hand, the chiral condensate is not a conserved quantity, and hence 
we cannot use the equation of continuity. However, this does not mean that there does 
not exist the reversible current. The typical example is the critical dynamics of liquid 
helium[TJ El HI El E]- In this case, we should solve the coupled equation of the order 
parameter, the entropy density and the superfluid velocity. The order parameter and 
the entropy density are non-conserved quantities, but the phenomenological equations 
contain the reversible currents, one of which is proportional to the superfluid velocity. 
Here, we derive the reversible currents coupled to the chiral condensate by using the 
projection operator method[TTl EE1 113 EES 113 EE3 EE3 Hi • 

In the projection operator method, the equation of motion of an operator A(x) is 
given by 

^-A(x) = / d 3 xW(zLA(x), AV)) ■ (A(x'), ^(x"))" 1 ■ A(x", t) 
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+ 



jf* ds J d 3 xW(/(x, 8 ), / f (x'> 0)) ■ (A(x'), ^(x"))- 1 • A(x", t - s) 



(F, G) = / o " f Tr I [e-^e-^Fe^G] , (26) 



+ /(x,t), (25) 

where the Liouville operator L is defined by L O = [H, O] with the Hamiltonian 
H §. The Hamiltonian introduced here is a microscopic Hamiltonian like the QCD 
Hamiltonian. However, we do not specify it because the Hamiltonian is irrelevant to 
define the reversible currents, as we will see later. Here, the first term on the r.h.s. is 
called the streaming term, the second is the memory term and the third is the noise 
term defined by /(x, t) = Qe lL ®( t ~ to ^iLA(iL,to). The Kubo's canonical correlation is 
defined by 

J 

where F and G are arbitrary operators, and Z = Tr[e~^ ]. 

The streaming term represents the contribution from reversible currents. Thus, in 
order to find the reversible currents coupled to the chiral condensate, we should find an 
operator O that satisfies the following condition, 

(iL*(x), O(x')) = (|<7(x), O(x')) = Tr[p[a(x), O(x')]] ? 0. (27) 

First, we assume O = vbTvb. Then, T = 7* and 7°7* can satisfy this condition, because 
[cr(x), O(x')] 7^ (This condition corresponds to the non-vanishing Poisson bracket in 
PP). However, the net baryon number current that is given by V = 7* is the reversible 
current for the net baryon number density. As a matter of fact, it is easy to calculate 
the expectation value in the free gas approximation. Then, one can see that Eq. (|27j) 
vanishes for T = 7*. It follows that the net baryon number current is not directly coupled 
to the chiral condensate. Thus, the reversible current of the chiral condensate can be 

J CT (x) = <Kx)7° T V(x) - $(x) 7 T Y(x))eff. (28) 

We are particularly interested in the coupling of the chiral condensate and the net 
baryon number current, because, as is discussed in [Hj, the net baryon number density 
is a promising gross variable with large time scale. However, if we put F = 7 , one 
can show that [cr(x), n(x')] = 0. Thus, the net baryon number density is not coupled 
to the chiral condensate directly. Now, we consider nonlinear couplings and assume 
O = (■07°"0) (ipT'i/j) . In order to satisfy the condition [cr(x), O(x')] 7^ 0, we should 
choose, for example, T' = 7°7*. 

Finally, the Langevin equation for the chiral condensate including the reversible 
current is given by 

<t(x) = - RV ■ J ff (x) - TV ■ (n(x)J CT (x)) - 7ss X ff (x) + f ff (x). (29) 

§ Of course, we can apply the projection operator method to classical field theory. Then, the Liouville 
operator is replaced by the Poisson bracket. However, for simplicity, we treat all gross variables as 
operators in the definition of the reversible currents, and we replace them with the corresponding 
classical variables at the last of the calculation. 
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It should be noted that this definition of the reversible coupling is only the necessary 
condition, because there is a possibility that the coefficients of reversible currents R and 
T vanish after calculating the expectation value in Eq. (|27j) . The coefficients depend 
on the choice of the microscopic Hamiltonian H and hence microscopic calculations 
are needed. However, in our phenomenological discussion, we leave them as unknown 
parameters. After all, to describe the critical dynamics near the chiral phase transition 
including the reversible couplings, we should solve the coupled equation of a, n, J n and 
J CT . 



5. Analysis based on the mode coupling theory 

In this section, we show that the reversible coupling derived above can affect the 
relaxation of the chiral condensate and change the dynamical universality class of QCD. 
First, we make the following ansatz for the chiral current; 

J CT (x) = a(x)v CT (x). (30) 
We further assume that the reversible flow of the chiral condensate are incompressible, 

V-v (T = 0. (31) 

Then, the Langevin equation for the chiral condensate is given by 

<t(x) = - i?K(x) • V)a(x) - T(v CT (x) ■ V)(n(x)a(x)) 

-7 ss X CT (x)+^(x). (32) 

The second term on r.h.s. gives rise to the coupling of the chiral condensate and the 
net baryon number density. However, this is a higher order nonlinear term and we drop 
it for simplicity. 

To estimate the effect of the nonlinear reversible coupling, we apply the technique 
developed in the mode coupling theory, where the nonlinear effect is calculated 
perturbativelyjH El EJ ||- The nonperturbative solution is calculated by setting R = 0, 

(T°(q, t) = e- iW2 V(q) + C dse- luJ2(t - s) £ a {q, s). (33) 

Jo 

We further introduce the time correlation function of the chiral condensate for t > t' > 0, 

G q (t - t') = (a(q, t)a(-q, t'))/ Xq , (34) 

where Xq = ( (T (q) cr ( — q))i an d ( ) denotes the expectation value for an initial state 
including the weight function of the noise. When we ignore the nonlinear term, the 
correlation function is 

Gj(f - = (trWt-q^OVXq = e"^-''). (35) 
With the help of the nonperturbative solution, the solution of Eq. (f3~2"|) is given by 

a(q, t) = a°(q, t) - R f tfci £ e^'-^k • v ff (q - k, s))<r(k, s), (36) 

Jo V k 

|| The mode coupling theory in the chiral phase transition is discussed also in ^1 
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Then, the time correlation function is developed by the following equation, 

jG^t) = -iu 2 G\{t) - -^ £<(ik ■ ~ k, t))a(k, t)a(- q )). (37) 

Substituting Eq. (f3l)J) into the second term on the r.h.s. of Eq. (fHTj) and expanding 
in powers of the R up to second order, we have 

jG^t) = - i^G^t) -^E((*.v ff (q-k, t))a°(k, t)o-(-q)) 



R 2 

+ 

k,k 



£ f dse-^- s \{^ ■ v a (q - k, t))(ik' • v CT (k - k', s))a°(k', s)a(-q)). 



We further apply the decoupling approximation as follows; 

((ik.v <T (q-k,t))<T°(k,t)(7(-q)> 

« ((iq- v CT (0,t)))(a (q,t)a(-q))^, (39) 
((ik • v CT (q - k, *))(ik' • v CT (k - k', s))a°(k', s)a(-q)) 

w <(ik • v ff (q - k, t))(iq • v CT (k - q, a))) (a°(q, £)a(-q))5g q . (40) 

Then, we obtain 

^G q (t) = - iu 2 G\{t) - |((iq ■ v CT (0,t)))G q (t) 

+ f Q ds^- 2 £ - 8)({ik ■ v CT (q - k, t))(iq • v CT (k - q, S ))>Gj(a). 

(41) 

This equation includes the reversible coupling up to second order. To resum the 
second order correction, we replace one G q (t) in each term with the full time correlation 
function G q (i), 

jG^t) = - tLU 2 G n (t) - 5((iq • v CT (0, *)))G q (t) + i? 2 jf* dsE^t, S )G q (s), 

(42) 

where we introduced the memory function H q (i, s), 

s ) = £<(* k • v -(q - k ' *))(<q • v -( k - q> *))>g,(* - *)■ (43) 

k 

This resummation corresponds to the ring-diagram approximation jH EJ. If we can 
employ the Markov approximation, that is, the memory function relaxes rapidly and 
vanishes at late time, the time dependence of the memory function is approximately 
replaced by the delta function, 

S q (t,s)«2H q 5(t-s), (44) 

with 

lim / dsEJt,s). (45) 

t— >oo Jo 
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We further assume a simple diffusion equation for the dynamics of the velocity field 
and ignore the reversible coupling [3 J, 

(|-^V 2 ) V(T (x,t) = 0. (46) 

Then, the correlation function of the velocity field is given by 

£<(ik • v CT (q - k, t))(iq ■ <(q - k, s))> = - £ 5>„(k, f) • v^k, s )> 
k J k 

= -^E^ k2( ^ s) (v CT (k)-v;(k)). 
6 k 

(47) 

Here, the initial correlation of the velocity is assumed to be Gaussian. Thus, we have 



q 2 ^ (v f (k)-v;(k)) _ (4f 



q W 2 Y -z/k 2 - icu 2 
Finally, we obtain the renormalized damping rate of the chiral condensate, 

Because of the critical slowing down [T5 | H6 | II 7j . iu)2 becomes smaller and vanishes as the 
temperature approaches the critical temperature ^f. If the velocity field also becomes 
slower as the temperature is lower to the critical temperature, and hence the diffusion 
constant v vanishes at the critical point as is the case with diffusion constants in spin 
diffusion processes and binary fluid mixtures^, the third term on the r.h.s. becomes 
larger than the first term, where the effect of the irreversible coupling is included. This 
is the reason why the reversible coupling can be more important than the irreversible 
coupling and the latter is ignored in various cases + . Then, the third term plays an 
important role in the critical dynamics of the chiral phase transition. This suggests that 
there is the possibility that the critical dynamics is changed by introducing reversible 
currents. 



6. Summary and concluding discussions 

We have discussed the phenomenological equations to describe the critical dynamics of 
the QCD phase transition. 

First of all, we derived the Langevin equations without the reversible couplings. 
We adopted the time- dependent Ginzburg-Landau equation for the time evolution of 
the chiral condensate and the Cahn-Hilliard equation for that of the net baryon number 
density. However, because of the Curie principle, the off-diagonal Onsager coefficients 
vanish. 

% This disappearance occurs in the chiral limit. When we have finite current quark masses, iuj-i can 
still be small but finite even at the critical temperature. 

+ Of course, we cannot say anything about the nonlinear irreversible coupling in this discussion. 



Phenomenological approach to the critical dynamics of the QCD phase transition revisited 11 

Next, we discussed the reversible coupling in the QCD critical dynamics. The chiral 
condensate is coupled not only to the net baryon number density but also to the chiral 
current through the reversible coupling. The reversible coupling plays an important role 
for the relaxation of the chiral condensate, because of the critical slowing down and the 
decreasing diffusion constant at the critical points. Then, the effect of the irreversible 
coupling is small [TJ El IH El E]- 

It follows that the dynamical universality class can be changed from model H 
discussed by Son and Stephanov. They claimed that, if we can ignore the coupling 
to the energy density and the momentum density, the universality class belongs to 
model B, and if we take them into account, it is changed to model H. 

In this study, we found that their discussion can be changed as follows: 1) Son 
and Stephanov concluded that the linear combination of the chiral condensate and the 
net baryon number density plays the role of an order parameter as in the case of the 
liquid-gas phase transition. However, if we take the possibility of the critical slowing 
down into account, the appearance of the reduction of the gross variables is not clear. 
2) The chiral condensate can be coupled to reversible currents, i.e. the chiral current 
and this coupling can change the dynamical universality class. 

If we take the reversible coupling of the chiral current into account, we must solve 
the coupled equation for the chiral condensate and the chiral current. The chiral 
condensate is not conserved and the chiral current is assumed to obey the diffusion 
equation. This situation is similar to the isotropic antiferromagnet, where the dynamical 
universality class belongs to model G. In the isotropic antiferromagnet, the order 
parameter is the three-component non-conserved vector density. The order parameter is 
coupled to a three-component conserved vector density that obeys a diffusion equation. 
In the chiral phase transition, the order parameter has four component in general, i.e., 
a and 7r\ Therefore, we can conclude that the dynamical universality class is similar to 
model G, noting, however, that the dimension of the order parameter is not three but 
four. 

In this calculation, we ignored the reversible coupling to the net baryon number 
density, the net baryon number density current and so on *. When we take them into 
account, the dynamical universality class can be changed again. Now, we have five 
gross variables that possibly couple to the chiral condensate; the chiral current, the net 
baryon number density, the net baryon number current, the energy density, and the 
momentum density. The several gross variables may not be independent as is the case 
with the liquid-gas phase transition. To discuss the possibility, we need a microscopic 
calculation as is done in ^Hl EH EE] • 

The criteria for choosing gross variables to describe the QCD critical dynamics are 
discussed in [2J EH] and the criteria do not allow reversible currents as gross variables. 
However, there is no general criteria to choose the complete set of gross variables jU El Ej- 

* The importance of conjugate fields as the candidates of gross variables near the chiral phase transition 
is discussed in X9 ■ Furthermore, it is known that if we ignore the conjugate held, it sometimes happen 
that we cannot separate slow and fast time scales [TCI II 7) . 
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In this paper, we have assumed that the chiral current is a gross variable. Of course, to 
see whether the reversible current is slow enough to be interpreted as a gross variable, we 
should calculate the time correlation function of the reversible current that characterizes 
the time scale of the reversible current [T7] . However, it will be natural to expect that the 
coarse-grained equation of the chiral condensate has the contributions of the reversible 
currents, because the transfer of the fluctuations of the chiral condensate should be 
induced by the macroscopic flows of quarks and gluons, even when the flows do not 
produce entropy. As a matter of fact, the fluctuations vary even in thermal equilibrium. 

We derived the Langevin equation to describe the QCD critical dynamics based on 
nonequilibrium thermodynamics, where the dissipation process is described by diffusion. 
However, as is pointed out in many papers, normal diffusion equations do not obey 
causality; their propagation speed is infinite and beyond the speed of light. To avoid 
this difficulty, extended irreversible thermodynamics (EIT) has been proposed where 
diffusion processes have finite relaxation time [2TJ 122 E31 EH • Thus, to discuss the 
QCD critical dynamics, it may be better to use EIT As a matter of fact, the microscopic 
calculation of the chiral condensate shows dissipation exhibiting oscillation and seems to 
support extended thermodynamics (THIIIZI HE]- Then, the critical dynamics in relativistic 
systems can be different from the non-relativistic cases ji|27j. 
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